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Abstract 

Main result of this paper is to derive the exact analytical expressions of 
information and covariance matrices for multivariate Burr III and logistic 
distributions. These distributions arise as tractable parametric models in 
price and income distributions, reliability, economics, populations growth and 
survival data. We showed that all the calculations can be obtained from one 
main moment multi dimensional integral whose expression is obtained through 
some particular change of variables. Indeed, we consider that this calculus 
technique for improper integral has its own importance in applied probability 
calculus. 
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1. Introduction 

In this paper the exact form of Fisher information matrices for multivariate Burr III 
and logistic distributions is determined. It is well-known that the information matrix 
is a valuable tool for derivation of covariance matrix in the asymptotic distribution 
of maximum likelihood estimations (MLE). In the univariate case for Pareto (IV) 
and Burr XII distributions, the Fisher information matrix expressions are given by 
Brazauskas |3] and Watkins 12 . As discussed in Serfling [TT], under suitable regularity 
conditions, the determinant of the asymptotic covariance matrix of (MLE) reaches an 
optimal lower bound for the volume of the spread ellipsoid of joint estimators. 

The univariate logistic distribution has been studied rather extensively and, in fact, 
many of its developments through the years were motivated to the normal distribution 
(see for example the handbook of Balakrishnan [3]). However, work on multivariate 
logistic distribution has been rather skimpy compared to the voluminous work that 
has been carried out on bivariate and multivariate normal distributions (Gumbel 8 , 
Arnold Johnson, Kotz and Balakrishnan J5] and Malik and Abraham ^H]). For 
a broad discussion of logistic models and diverse applications see Malik and Abra- 
ham [TH] , 

Burr III and Logistic distributions also arise as tractable parametric models in 
the context of actuarial science, reliability, economics, price and income distributions 
(Dagum [7], Burr [o] and Burr |oj). 

This paper is organized as follows: Multivariate Burr III and logistic distribution 
are introduced and presented in sectional Elements of the information and covariance 
matrix for multivariate Burr III distribution is derived in section |2] Elements of the 
information and covariance matrix for multivariate logistic distribution is derived in 
section 0] Conculusion is presented in section [5] Derivation of first and second deriva- 
tives of the log— density function of multivariate Burr III distribution and calculation 
of its main moment integral are given in Appendices A and B. Derivation of first and 
second derivatives of the log density of multivariate logistic distribution and calculation 
of its main moment integral are given in Appendices C and D. 
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2. Multivariate Burr III and logistic distributions 

The density function of the Burr III distribution is 

where — oo < fi < +00 is the location parameter, 9 > is the scale parameter, c > is 
the shape parameter and a > is the shape parameter which characterizes the tail of 
the distribution. 

The n-dimensional Burr III distribution is 

(\ — (ct+n) 

where x = [xi,--- ,x n ], Xi > fa, Cj > 0, —00 < fa < +00 , a > 0, 0j > for 
i = 1, • • • , n. One of the main properties of this distribution is that, the joint density 
of any subset of the components of a multivariate Burr III random vector is again of 
the form [HI- 

The density of the logistic distribution is 

fx{x) = - e e-^(l + e-my {a+X \ x>», (3) 

where —00 < fi < +00 is the location parameter, 9 > is the scale parameter and 
a > is the shape parameter. 

The n-dimensional logistic distribution is 




e~y^^J, (4) 

I-! 

i = l 



where x — [xi, ■ ■ ■ , x n ], Xi > fa, a > 0, —00 < fa < +00 and 0, > for i =!,■■■ ,n. 
The joint density of any subset of the components of a multivariate logistic random 
vector is again of the form (@J 9 . 



3. Information Matrix for Multivariate Burr III 



Suppose X is a random vector with the probability density function /e(.) where 
= (6*i, 9-2, 9k)- The information matrix 1(0) is the K x K matrix with elements 

rd 2 lnf e (Xy 



-E 



d9,d9, 



i,3 



1. 



■ K. 



(5) 
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For the multivariate Burr III, we have 9 = (/Ui, • • • , [i n ,9i, • ■• , n , c\, ■ ■ ■ ,c n , a). In 
order to make the multivariate Burr III distribution a regular family (in terms of 
maximum likelihood estimation), we assume that vector /j, is known and, without loss 
of generality, equal to 0. In this case information matrix is (2n + 1) x (2n + 1). Thus, 
further treatment is based on the following multivariate density function 




-(a+n) 



f[ {a + i - 1)Ci ( X A iCi+1 \ ^>0- (6) 



i=i 



The log-density function is: 

n 

In f n (x) = + i - 1) - + In d\ - (<* + 1) In 



i=l 



(a + n) In 1 + £ 



(7) 



Since the information matrix 1(0) is symmetric, it is enough to find elements 7^ (O), 
where 1 < i < j < 2n + 1. The first and second partial derivatives of the above 
expression are given in the Appendix A. In order to determine the information matrix 
and score functions, we need to find expressions of the generic terms such as 



E 




and 



E 



0, 



and evaluation of the required orders partial derivatives of the last expectation at the 
required points. 

3.1. Main strategy to obtain expressions of the expectations 

Derivation of these expressions are based on the following strategy: first, we derive 
an analytical expression for the following integral 



E 



+oo 



n ? /«(*) d *> 



(8) 



and then, we show that all the other expressions can easily be found from it. We 
consider this derivation as one of the main contributions of this work. This derivation 
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is given in the Appendix B. The result is the following: 

'Xi 



E 



n 



OC Z' + OC 

ii Jo 



n 

i=l 



/ n (jc)da; 



r(a) 



where T is the usual Gamma function and 



r(q + £IU*)II?=ir(i-%) ^ r, r, 

— = > — < a, — < 1, 



i=l 



(9) 



a 



E 



i=i a ) 
-= ^, 1 



dr k dn 



< Lk < n, 



*W( z ) = ^L/' r/ W 



z > 0, 



dz" V r ( z ) 
9(m+n) flnnXfA = fl&("*+n)f \ > o 



9r' 

and n, m £ N. 

Specifically, we use digamma trigamma $f'(z) and \E , rirfc (z) func- 

tions (Abramowitz pQ and Brazauskas 0]). To confirm the regularity of ln/ n (x) and 
evaluation of the expected Fisher information matrix, we take expectations of the first 
and second order partial derivatives of J7J). All the other expressions can be derived 
from this main result. Derivative with respect to a, from the both sides of the relation 



+ OC 



fn(x) dx, 



leads to 



E 



n 1 

y — - — • 

f-' a + i - 1 



(10) 



(11) 



V W*l r (« + 2 + ^) r (-S + !) r (-S + 1) 



From relation ©, for a pair of (I, k) we have 

f(ri,r k ) = E 
From relation 112|) . at r k =0we obtain 

'XtY'] r(n+a)r(-n+i, 



r( tt ) 



(12) 



C(n) = E 



0/ 



r(a) 



(13) 
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Evaluating this expectation at r; = — cj, r; = — 2q and the relation l|12[) at (n = 
-Q,r fe = -Cfc), we obtain 



E 



E 



X, 



X, 



-2c, 



E 



Xi\- Cl [X k ^ Ck 



a-T 



(a-l)(a-2)' 
1 

(a-l)(a-2)- 



Evaluating the required orders partial derivatives of 113( 1 and (|12fl at the required 
points, we have 




[*(<*)- r'(i)] 

Q 

[*(<*- 1)-T'(2)] 
cjr(a-l) 

[* 2 (a - 1) + *'(o - 1) - 2*(a - l)r'(2) + r"(2)] 
Q 2 (a-1) ' 

[2* 2 (a - 2) + 2*'(a - 2) - 2#(a - 2)F(3) + r"(3)] 
q 2 (a-l)(a-2) 

[2^(a - 2) - r(3)] 
d(a-l)(a-2) ' 

[*Ma- 2) -T'(2)] 
Ck(a-l)(a-2) ' 

[-r(2)(W r ,(a-2) + ^ rfc (a-2))] 
Qc fc (a - l)(a - 2) 

[(E(2)) 2 + ^ ;rt (a-2)] 
Qc*(a — l)(a — 2) 



From these equations with a replaced by (a + 1) and (a + 2) in (jSJ, we can show 
that 
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E 



E 



(i + 

_(i+E; =1 (t)- c 2 




,(i+E; =1 (f)- 
(£)- 2c; Kf) 



E 




E 



_ (i+E; =1 (f)-^) 

(f)--(f)--ln(f) 



(i+E; =1 (tH" 



-E 



a + n 



a+l 



1 



a + n 



a(a + 1) 



(a + + n + 1) 
2 



E a +2 



X, 



-2c, 



(a + + n + 1) ' 

[* 2 (a) + tf'(a) - 2*(a)r'(2) + r"(2)] 
Q 2 (a + n) 



[*(<») -T(2)] 
c;(a + n) 

[2^(a)-r(3)] 
Q(a + n)(a + n + 1) ' 

[2* 2 (a) + r"(3) - 2r'(3)*(a) + 2*'(a)] 
c 2 (a + n)(a + n + 1) 

[*,-»- r(2)] 

c k (a + n)(a + n + 1)' 

[tt r , r >) - Hg) QMa) + *Ma)) + (r(2)) 2 
Qc k (a + n)(a + n + 1) 



3.2. Expectations of the score functions 

The expectations of the first three derivations of the first order follow immediately 
from the corresponding results for their three corresponding parameters and we obtain: 



E 



E 



dlnf n (X) 
da 



£ 

»=i 



a + i — 1 



E 



\ din f n (xy 


Cl 


( (a + n)ci \ 


[ de l 


~ 0i 


I Oi J 



ln ( 1+ E(^) 

j=i j 

X±\-ci 



J \ -Ci 



E 



(If) 

(i+e; =1 (^)- c 



= 0, 
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E 



dln/»(X) 
dci 



= E 



In 



X, 



+ (a + n) E 



i + E- =1 (^)- 



= o. 



3.3. The expected Fisher information matrix 



Main strategy is again based on the integral which is presented in the Appendix 
B. After some tedious algebric simplifications, the following expressions can be ob- 
tained 



Ix{a) = 

Ix{6i,a) = 

Ix(ci,a) = 

lx(0i) = 

Ix(ci) = 

Ix(0i,6k) = 

Ix(ci,c k ) = 

Ix(0hC k ) = 

IxtfhcA = 



E 



i 



i (a + i — 1)' 

Cl 



9i (a + n)' 
1 

ci (a + n) 
cf(a + n — 1) 
(9 2 (a + n+ 1)' 
1 



I = 1, ■ •• ,n, 



[*(a) - T'(2)] , l = l,...,n, 



1 + T"(2) - 2*(a)r'(2) + * 2 (a) + *'(a)] 
r'(3)*(a) + ^ 2 {a) + *'(a) 



cf (a + n + 

CjCfc 

0/0 fc (a + n + 1)' 



1 = 1, 



-r"(3) 



i = i. 



(r'(2)) 2 - T'(2) (* rfc (a) + * ri (a)) + (a) 



QCfe (a + n + 1) 



[-r(2)+^(a)] 
9ic k (a + n + l) ' 
[F(2) -*(«)] [2^(a)-r(3)] 
0, 0/ (a + n + l) • 



/ = !,••• ,71. 



. n, 



Thus the information matrix, / B U1T ni(®)' ^ or ^ ne multivariate Burr 111(0, 9, c, a) 
distribution is 



1 Burr IIl( e ) 



I{6 u 6 k ) I{9 u c k ) I(6 h a) 
I(0i,c k ) I(ci,c k ) I(ci,a) 
I(6 u a) I{c u a) 1(a) 



(14) 
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3.4. Covariance matrix for multivariate Burr III 

Since the joint density of any subset of the components of a multivariate Burr III 
random vector is again of the form (0) we can calculate the expectation 



E 



v \ mi / y 



oo poo / \ m i / \ m k 

' xi- m\ / x k - fi k 



0! 



0k 



fx u x k (xi,x k )dxi dx k 



T( a + 2ii + m.)T(i - m)r(l - no.) 

v ci e fc I V Cj - 1 V c k > 

r(a) 

m; mi. m; mi 

1 ->0, 1 ->0, a+ — + — >0. 

Q c fe c; c fc 



(15) 



Evaluating this expectation at (m/ = 1, m k = 0), (to/ = 0, = 1) (m; = 1, 
rnfe = 1) and (nij = m, rrife = 0) we obtain 



E[Xi] = 



E [X k ] 
E [X,X fe ] 

Var{X ; } 
Cov [X h X k ] 



T(a) ci ci 

/ifc + =^r[T(a+i)r(l-i)], 
r(aj Cfe c fe 



= ^ifcE [X;] + /i;E [Xk] - Hl^k 

+^[r( a + - + -)r(i-!)r(i-I)], 

r(aj Cfe Cfc Cfc q 

= =f T p(a + -r(i--) ) 

I (a) c; q 

r(a + -)r(i - -)i» - r 2 (i - -)r 2 (a + -) 

Cl Ci Cl Cl 



T 2 (a) 

01 k 



1 > 0, 

Cl 



-r(i - -)r(i - l)[r(a + - + -)i» - r(a + -)T{a + i; 

r 2 (a) Q Cfc. q Cfc Cfc ci 



1 < I < k < n, 



fc = 2,--- 



4. Information Matrix for Logistic distribution 



For the multivariate logistic distribution, we have = (fii, ■ ■ ■ , fj, n , $x, ■ ■ ■ ,0 n ,a). 
In order to make the multivariate logistic distribution a regular family (in terms of 
maximum likelihood estimation), we assume that vector /j, is known and, without loss 
of generality, equal to 0. In this case information matrix is (n + 1) x (n + 1). 
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Thus, further treatment is based on the following multivariate density function 

/ \ -(a+n) 

/„(.)- [l +£.-(*)] f l <°±£H e -W. (10) 



Thus, the log-density function is: 



lnf n (x) = J2 N« + i ~ 1) - Inft] - ^(|) " (« + «) ln ( 1 + E e ~ ( ^ ) ' ( 17 ) 



Since the information matrix is symmetric it is enough to find elements 
where l<i<j<n+l. The first and second partial derivatives of the above 
expression are given in the Appendix C. Looking at these expressions, we see that to 
determine the expression of the information matrix and score functions, we need to 
find the following expectations 



E 



and E 



< v "' I I e V e k 



and evaluation of the required orders partial derivatives of the last expectation at the 
required points. 

4.1. Main strategy to obtain expressions of the expectations 

Derivation of these expressions are based on the following strategy: first, we derive 
an analytical expression for the following integral 



E 



n 



p + OO p + OO " / / x .\ \ Ti 

l - l n^ M ) aw* 



(18) 



and then, we show that all the other expressions can be found from this easily. This 
derivation is given in the Appendix D. The result is the following: 



E 



n(e-w) = / -J II 

<=1 v ' J ° J ° i=l 



f n (x)dx 



r(a) 

^V;<a, ri>-l. (19) 



i=l 
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Taking of derivative with respect to a, from the both sides of the relation 

r+oo 



leads us to 



E 



1 = 



ln[l + £, 

i=i 



f n (x)dx, 







n 

y — - — • 

^— ' a + i — 1 

i— 1 



From relation H19|) . for a pair of (I, k) we have 



<p(ri,r k ) = E 



e v 



r(a - n - r fe )r(r fe + i)r(n + i) 



From relation i|22|) . at r& = we obtain 

' r(a - ri)r(r, + 1) 



E 



T(a) 



and evaluating this expectation at r; = 1, we obtain 

1 



E 



e v "i 



a - 1 



(20) 



(21) 



(22) 



(23) 



(24) 



Differentiating first and second order of 112MI) with respect to n and replacing for n = 0, 
H = 1 and r; = 2, we obtain the following relations: 



E 



E 



AO 



E 



ft 



E 

From relation (|22|) . 
9 



dndr k 



*(a- 1) -T'(2) 



(a-1) 

* 2 (a - 1) - 2F(2)¥(a - 1) + *'(a - 1) + r"(2) 



(a-1) 

[r"(3) - 2r'(3)*(a - 2) + 2* 2 (a - 2) + 2#'(a - 2)] 



(a-1) (a -2) 



<p(ri = l,r k = 1) 



E 



Xl \ ( At 



g^(2) (r'(2)-$ ft (a-2)-$ r ,(a-2))] 
(a-l)(a-2) 



JVja_-2)_ 
(a-l)(a-2)' 

With a replaced by (a + 1) and (a + 2) in (JTBJ, we obtain 



(25) 
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E 



E 



E 



E 



(i + E 




(i + E 

(*: 




(i + E 

(*) 


~.<*>) 



E 



fL) e -(lf) e -(^) 



a + n 



a 



a + n 



E a +i 



E q +i 



1 



a + n 



X, 



n e -(9f) 



a + n 



^ 2 (a) + gg(a) + E"(2) - 2E'(2)^(a) 
a + n 

[2^ 2 (a) - 2r'(3)tf (a) + 2*'(a) + r"(3)] 
(a + n)(a + n + 1) 

r'(2)[E'(2) - ^ rk (a) - VE r , (a)} + gy, (a) 
(a + u)(a + n + 1) 



i+e;=i^~ 

4.2. Expectations of the score functions 

The expectations of the first two derivations of the first order follow immediately 
from the corresponding results for their two corresponding parameters and we obtain: 



E 



dlnf n (X) 
da 



E 



a + i — 1 



— E 



ln(l + £e 

3=1 



= o, 



E 



dhif n (X) 
09i 



1 + \ [*(«) - T'(l)] - i [^(a) T'(2)] = 0. 
v\ til t>i 



4.3. The expected Fisher information matrix 

Main strategy is again based on the integral (|18|) which is presented in the Appendix 
D. Again after some tedious algebric simplifications, the following expressions can be 
obtained 
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Ix{ot) 
Ix{Oua) 



E 



i 



—[{a + i-l) 
1 



Oi (a + n) 



[*(«)- T'(2)], 1 = 1 



1 } 3 



' l = ^(aVn + [*'(<») ~ 2r'(2)fr(a) + r"(2) + *'(a)] 



F(2)-*(a) 



/ = !,-•■ ,71, 



Ix(6i,0k) 



f(a + n+l) 
r'(2)[r'(2) - * r » - * n (a)] + * r;rfe (a) 



fc ^ I. 



6i6 k (a + n + l) 

Thus the information matrix, for the multivariate logistic (0, 8, a) distri 

bution is 

I(6 u 6 k ) I(6 h a) 
W,a) 1(a) 



'ML(e) = 



(26) 



4.4. Covariance matrix for multivariate Logistic 



Since the joint density of any subset of the components of a multivariate logistic 
random vector is again multivariate logistic distribution^), we can use the relation 
with g r ®g rk ( p( r i — 0) r k — 0) and obtain 



6,6. 



(T'(l)f - r'(l) (9 rk (a) - 9 n (a)) + * nrk (a) 



, k^l, 



E [XiXk] 

ELY;] = 6i [9(a) - r'(l)] , Z=l,...,n, 

E = fc [tf (a) - r'(l)] , fc = 1, • • • , n. 

From second order derivative of relation l|22|l . i.e., -^<p(ri = 0,r k = 0) we have 

E[Xf] = #[(r''(l))-2r^l)*(a)+* 2 (a) + *'(«)], l = l,---,n, 

Cov[X h X k ] = 9 l k [-r'(l)(9 rk (a) + 9 ri (a)~29( a )) + 9 nrk (a)-9 2 (a)] , fc ^ Z, 

Var{X,} = 6f[T"(l)-(T'(l)) 2 +9'(a)], l = l,.--,n. 



5. Conculusion 



In this paper we obtained the exact forms of Fisher information and covariance 
matrices for multivariate Burr III and multivariate logistic distributions. We showed 
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that in both distributions, all of the expectations can be obtained from two main 
moment multi dimensional integrals which have been considered and whose expression 
is obtained through some particular change of variables. A short method of obtaining 
some of the expectations as a function of a is used. To confirm the regularity of the 
multivariate densities, we showed that the expectations of the score functions are equal 



In this Appendix, we give in detail, the expressions for the first and second 
derivatives of ln/„(a;) , where, f n (x) is the multivariate Burr III density function 
which are needed for obtaining the expression of the information matrix: 



to 0. 



Appendix A. Expressions of the derivatives 




d\nf n (x) 
ddi 



ci (a + n)ci 




d\nf n (x) 
dci 




m- ci in (ft) 



l = !,-■■ ,n, 



(i + E? =1 (ff)- c 0' 



& 2 lnf n (x) 
dcida 



d 2 lnf n {x) 
ddida 





I = !,-■■ ,71, 
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d 2 lnf n (x) 
da 2 

d 2 lnf n (x) 
36? 



d 2 lnf n (x) 
dd 2 



n 1 

2—1 



Q (g + n)(l-Q) Q (ff)- c ' 



i+E; =1 (if)- c O 



+ 



(a + n)c 2 



( xn-2c, 



(i + E3Li(^)- 



! (ff)^ In 2 [f\ 
—~ - (a + n) t — - \ + (a + n) 

,f (i+E; =1 (it)^) 



1, • • • ,n, 

i X i)- 2Cl 1^ (ft) 



d 2 \nf n (x) = ( (a + n)c lCk \ (ff)- c '(ff)' 

m ' (i+e;-i(Ij) 



a 2 ln/ n (aQ 
dc k d8i 



dckdci 
d 2 \nf n (x) 



-, k ^ I, 



(a + n)cA v e J 



0),. 



= (a + n) 



(i+E; =1 (it)- c 

(ftr Cl (t)--ln(t)ln(ff) 



+ 



+e; =1 (^)- c ^ 
(i+E; =1 (t)^ 



1 / a + n 



0i J (i+E;=i(it)^) 

(a + n)cA (ff)- 2C!ln (lr 



k^l, 



0i 



(i + E^df)-^)' 



Appendix B. Expression of the main integral 

This Appendix gives one of the main results of this paper which is the derivation of 
the expression of the following integral 



E 



.i=i 



f+oo p+oc n , \Ti 



(i) 
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where, f n (x) is the multivariate Burr III density function ij^J. This derivation is done 
in the following steps: 

First consider the following one dimensional integral: 



+ OC 



aci I x\ 



"(ci + l) 



— (a+n) 



3=1 



+ OC 



ac\ I x\ 



r i \ -( c i + !) 



— (a+n) 



J=2 



(*)' 



-(a+n) 



dxi. 



Note that, goings from first line to second line is just a factorizing and rewriting the 
last term of the integrend. After many reflections on the links between Burr families 
and Gamma and Beta functions, we found that the following change of variable 

1 + ^ —1=7^-, 0<f<l, 



simplifies this integral and guides us to the following result 



(2) 



-(a+n)— n+1 



r(a + n) 



3=2 



(3) 



Then we consider the following similar expression: 

"+ 00 c 2 a{a + l)r(l - %)T(a + n - 1 + ^) 



C 2 = 



•T2 



6 2 T(a + n) 

"(C2+1) 



i=2 



-(a+n)-^-+l 



da; 2 , 



c 2 a(a + l)r(l - a)r( a + „ - 1 + n.) 



6< 2 r(a + n) 

-(cj+l) 



-(a+n)-|- + l 



_( a+n )_rt+i 



dx 2 , 
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we obtain: 



Co = 



i)r(i 



£ 

3=3 



ir)r(i 



-)r(a- 

r(a + n) 

_( a+n )_rL_ 



ci 



-) 



L +2 



(5) 



Continuing this method, finally, we obtain the general expression: 



C n — E 



n 



r(a) 



£ 



< a, 



>0. 



(6) 



We may note that to simplify the lecture of the paper we did not give all the details 
of these calculations. 



Appendix C. Expressions of the derivatives 

In this Appendix, we give in detail the expressions for the first and second derivatives 
of ln/ n (a;), where, f n (x) is the multivariate logistic density function (|16|1 . which are 
needed for obtaining the expression of the information matrix: 
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dhif n (x) 
da 

d\nf n (x) 
86i 

d 2 \nf n (x) 
da 2 

d 2 lnf n (x) 
86 id a 

d 2 \nf n (x) 



4^ a + 1 - 1 \ f-f 



(a + n) 













n 1 

l — l V 7 



z = i 



2 ' 7^ ') 



d 2 ln/ n (a;) 1 2a; ; 2 (a + n) 



(*) 



e -(^) 



of of 



(a + n) 



1 + E" =1 e^ 



(*)' 



(a + n) 



i + E*U 



/ = i, 



,n, 



(t) ; 



e e i 



Appendix D. Expression of the main integral 

This Appendix gives the second main result of this paper which is the derivation of 
the expression of the following integral 



E 



n 



+oo r+°° n 



n 



f n (x)dx, 



(1) 



where, f n (x) is the multivariate logistic density function This derivation is done 

in the following steps: 
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First consider the following one dimensional integral: 

-(a+n) 

dx n 



e 

3=1 



-(a+n) 



dx n . (2) 



Note that, goings from first line to second line is just a factorizing and rewriting the 
last term of the integrend. After looking for the links between logistic function and 
Gamma and Beta functions, we found that the following change of variable 

simplifies this integral and guides us to the following result 

(\ — (a+n)+r n +l 
i+i.-^i . <4) 
i^i ) 

Then we consider the following similar expression: 

r _ r(r n + l)r(a + n-r n -l) [ +oc (a + n - 2) / - ( -^ r ' 
" T(a + n-l) J 

-(a+n)+r n + l 

n—1 



r(r„ + l)r(a + n - r n - 1) /"+ 00 (a + n - 2) / _ ( |»=i) 



JO 



r(a + n-l) y 6 n -i 

\ -(a+n)+r„ + l / «„_! \ -(a+n)+r„4 



and again using the following change of variable: 
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we obtain: 



C„_i — 



r(r w _i + l)r(r» + l)r(a + n - r n - r n _i - 1) 
T(a + n- 2) 

-(a+7l)+r n +r„_i + l 



n-2 



e j 



(6) 



Continuing this method, finally, we obtain the general expression: 



Ci = E 



n 

.i=i 



r(«-EL,r ( )iE. 1 r(i ±2 ) i £ n<a _ 



r(«) 



i=l 



We may note that to simplify the lecture of the paper we did not give all the details 
of these calculations. 
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